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FOREWORD 


This  report  documents  a  study  accomplished  during  the  period 
15‘ January  1970s  to  14  January  1971  by  the  Department  of  Aerospace  and, 
Mechanical  Engineering,  University  of  Notre  Dame,  Nofcre  Dame,  Indiana 
46556,  under  Contract  F08635-70-C -0033  with  the.  Air  Force  Armament 
Laboratory,  Air  Force  Systems  Command,  Eglin  Air  Force  Base, 

Florida  3254-2.  Program  monitor  for  the  Armament  Laboratory  was 
Dr.  Mark  0.  Glasgow  (dlgr). 


This  report  has  been  reviewed  and  is  approved. 


CHARLES  K.  ARPKE,  Lt  Colonel,  USAF 
Chief,  Weapons  Effects  Division 
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ABSTRACT 


Wind  tunnel  tests  were  conducted  to  study  the  effects  of  a  boattail 
on  the  dynamic  stability  of  the  three-caliber  Army-Navy  basic  spinner 
projectile.  Three -degrees -of-freedom  angular  data  was  obtained  in  vertical 
subsonic  wind  tunnel  tests  on  one  model  with  no  boattail  and  on  one  model 
with  a  7.5-degree  boattail.  Data  reduction  yielded  the  restoring,  damping, 
and  Magnus  moments  as  functions  of  time  and  nonlinear  functions  of  angle 
of  attack  for  each  model.  The  testing  technique  employed  a  low  friction 
jewel  bearing  model  mount  and  a  means  for  obtaining  two-axis  data. 
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SECTION  I 
INTRODUCTION 


The  design  of  small  arms  projectiles  has  improved  from  the 
classical  round  shot  to  the  elongated  cylinder  aft  of  a  pointed  nose.  A  need 
to  obtain  greater  range  and  improved  accuracy  has  prompted  these  con¬ 
figurational  changes  over  the  years. 

The  range  of  a  projectile  mainly  depends  upon  two  factors:  the  mass- 
to-frontal-area  ratio  and  drag.(^)  Elongation  of  a  sphere  into  a  cylindri¬ 
cal  form  greatly  decreases  the  frontal  area  and  also  allows  the  addition  of 
a  large  amount  of  mass  to  the  body.  Drag  is  also  decreased  by  virtue  of  the 
streamlined  shape. 

Still  less  drag  and ,  therefore,  greater  range  can  be  attained  by 
changing  the  projectile  to  the  boattail  configuration. w  However,  a  boattail 
can  be  the  source  of  dynamic  instability  in  a  projectile  and  the  increase  in 
drag  due  to  the  large  pitching  and  yawing  motion  associated  with  unstable 
flight  can  overcome  the  advantage.  Therefore,  it  is  imperative  to  deter¬ 
mine  what  aerodynamic  effects  are  caused  by  a  boattail  in  relation  to 
dynamic  stability.  This  report  presents  a  comparative  study  of  the  dynamic 
stability  characteristics  of  a  basic  research  configuration  with  and  without 
a  boattail. 
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SECTION  II 


AEROBALLISTIC  THEORY 


For  constrained  center  of  gravity  wind  tunnel  tests-,  the  normal 
force  contribution  is  zero.  The  aerodynamic  forces  and  moments  are 
shown  in  Figure  1. 


Z 


Figure  1.  Aerodynamic  Forces. and  Moments 
Linear  Theory 

If  linear  variations  with  angle  of  attack  are  assumed,  the  aero¬ 
dynamic  moments  can  be  written  as 

M+iN=-iMaa  +Mqq-iM^a  -PMpaof-M^8f  e ipt 

The  resulting  complex  differential  equation  of  motion  is 
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lol 


Mpa  (|al>  ~  MP<*0  +  MPa2 


(The  Nonlinear  Theory  used  assumes  an  approximate  solution  of  the  same 
form  as  the  Linear  Theory;  however,  the  stability  parameter  now  reflects 
the  nonlinearity  in  the  stability  coefficients.) 

The  solution  for  the  complex  angle  of  attack  is  given  by^ 


0  =  KNe  i<^>Nt  +  Kpe  l^t 

(14) 

0  =  0  +  i^ 

(15) 

Kn  =  KNq  .  KP  =  KPq  eXY 

(16) 

X*N,P  =  XN,P  +  XVN,P 

(17) 

xv  _  QSd2 
N,P“  ~W 


2 IV 


\2  c:ma  r  n  ,  P  >p  +  2K p Pf  N) 


=  QStf 


N,P"  2V 


dCpl^r  m<*2 

[Cmq  +  Cm  a]  [Cmo  +  Cma] 

i - - ±o  n+  m  p)  +  - - - 


(18) 
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[kn,p+kp,N](1±n,p)  +  kp,ntn,p(p>n  1  x) 


+  c  „  +  C~-  -  LnI? 


-  mPao  mP«2  *  N,P  — 


(19) 


P*x 


N.P  "  ~oT 


(1  +  _L  )t  +  0N  ,P 


(20) 


STABILITY  GOEF F.ICIE NTS 


To  obtain  aerodynamic  stability  coefficients  from  the  angular  data 
obtained  from  the  wind  tunnel  tests,  the  Aeroballistic  Theory  is  fitted  to  the 
^angular  orientation  data,  9,<P  .  This  is  done  through  the  use  of  the  "Wobble” 
computer  program.^)  This  program  fits ‘the  theory  to  short  segments  of 
the  data  in  overlapping  sections  so  that  the  stability  parameters,  Kjsj,p, 

*N  P  >  “N  P  »  are  determined  as  functions  of  time.  Since  the  Nonlinear 
The’ory  assumes  an  approximate  solution  of  the  same  form  as  the  Linear 
Theory,  the  "Wobble"  program  can  be  applied  to  both. 


Computation  of  Linear  Coefficients 

Using  Xm  p  and  p,  along  with  the  velocity,  dynamic  pressure, 


roll  rate,  and  physical  parameters  of  the  projectile,  the  aerodynamic 
stability  coefficients,  Cm  ,  C  +  Cm£  ,  and  C  , 
functions  of  time  from  the  following  equations: 

81 
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6Jp 

rrpd  & 


C  + 
mq 
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=  + 


>x 


mp  a  —  — 
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2V 


=  2IV 
QSd2 


+  ^  p) 


Pmq-Cm^Xl 
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(26) 


(27) 


(28) 


Computation  of  Nonlinear  Coefficients 


The  nonlinear  aerodynamic  stability  coefficients,  Cm^  (  a  ), 
Cmq  (  «  )  +  Cmi(  <?),  and  Cmpa  ( ~a  )  were  computed  as  polynomial 
functions  of  the  complex  angle  of  attack  as  follows: 


a>Nu>p= 
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o  + 
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QSd 


(29) 
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where 


9  9  2 

8Za  =  Kn  +  Kp  + 


Knwn 


-  Kd<u 


P"  p 


-  Ct>, 


(30) 


Using  a  least  squares  technique  to  fit  a  straight  line  to  versus  Ba 

yields  CmQ,^  as  the  intercept  and  CmQf2  as  the  slope. 

Correcting  X*j  and  X*  by  determining  ^  and  Xp  from  a 

logarithmic  technique  developed  in  Reference  6,  the  actual  damping  rates 
XN  and  Xp  are  fitted  simultaneously  with  a  least  squares  procedure  to 

yield  (Cmq+  Cmcj)0  >  (cmq  +  Cmo^2*  ^poo’  and  Cmpa2* 
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SECTION  IV 

'  THREE -DEGREES-OF-FREEDOM 
WIND  TUNNEIi  TESTING  TECHNIQUE 


A  unique  technique  was  developed  for  obtaining,  three -degrees -of - 
freedom  angular  data  for  the  Army-Navy  basic  spihner  in  a- subsonic  verti¬ 
cal  wind  tunnel.  The  test  model  was  mounted  on  a’ sting. with  a  needle¬ 
point  set  into  an  inverted  jewel  cup  bearing’ which  allows  roll,  pitch,  and 
yaw  with  a  minimum  of  friction'®' .  An  optical  reference  system  which 
introduces  absolutely  no  air  flow  interference  was  also  used. 

Motion  pictures  of  the  angular  motion  of  the  model  were  taken  with 
a  constant- frame  rate  camera- positioned  almost  vertically  above  the  model 
The  departure  of  the  line -of- sight  from  the  vertical  was  accounted  for  in 
the  data ‘reduction  process-.  The  data,  obtained  from  the  films  through  the 
use  of  an  optical  comparator,  was  converted  to  Euler  angles  and  fitted 
with  the  "Wobble"  fitting  program. 


Moment  of  Inertia 

The  axial  and  transverse  moments  of  inertia  were  found  with  the 
use  of  a  torsional  pendulum. 

The  axial  moment  of  inertia  was  theoretically  computed  for  a 
cylindrical  mass  of  known  density  and  dimensions.  The  mass  was  then 
suspended  from  a  stiff  wire  and  given  an  initial  angular  displacement. The 
period  of  the  oscillation  was  measured,  using  a  photo-diode.  Dividing  the 
period  of  oscillation  by  the  theoretical  moment  of  inertia  yields  the 
torsional  spring  constant. 

Knowing  the  torsional  spring  constant,  the  moment  of  inertia  for 
the  model  was  then  determined  by  suspending  the  model  from  the  mass  and 
again  measuring  the  period  of  oscillation  with  the  photo-diode.  The  total 
moment  of  inertia  minus  the  theoretical  moment  of  inertia  yields1  the 
moment  of  inertia  for  the  model.  This  was  done  to  find  both  the  axial  and 
transverse  moments  of  inertia. 


Data  Acquisition 

The  tests  were  performed  in  a  vertical  subsonic  wind-tunnel  which 
allows  the  use  of  a  low  friction  jewel  type  bearing. 
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The  jewel  bearing  is  positioned  securely  inside  the  model  on  its 
axis  of  symmetry  and  coincident  with  its  center  of  gravity.  The  model  is 
placed  on  a  needle  tipped  sting  of  sufficient  thickness  to  limit  vibrations 
which,  in  turn,  is  positioned  in  a  holder  mounted  inside  the  wind  tunnel. 
The  model  has  a  tiny  white  dot  painted  on  the  nose  for  reference. 

With  the  model  in  place,  a  high  speed  camera  which  operates  at  a 
constant  128  frames  per  second  is  positioned  outside  the  wind  tunnel  and 
above  the  model  (Figure  2).  (The  ideal  situation  would  be  to  have  the 
camera  looking  straight  down  at  the  nose  of  the  model;  however,  the  angle 
of  the  camera  with  respect  to  the  flow  is  taken  into  account  when  reducing 
the  data.)  The  movie  camera  is  then  focused  on  the  white  dot  on  the  nose 
of  the  model.  The  camera  is  equipped  with  a  lens  having  a  focal  length  of 
20  mm.  The  short  focal  length  is  needed  to  insure  a  large  depth  of  field 
since  the  nose  of  the  model  travels  slightly  above  and  below  the  plane  of 
focus  during  the  test. 


•-Camera 


*-*  Test 
Section 


Reference 

Source 


Figure  2.  Position  of  Test  Equipment. 


A  reference  system  is  employed  which  causes  no  flow  obstruction 
whatsoever.  To  obtain  the  true  motion  of  the  lest  model,  it  is  important 
that  nothing  be  placed  in  the  flow  so  as  to  alter  the  aerodynamics. 

A  black  box  with  two  pinholes  0.035  inch  in  diameter  and  with  a 
light  source  inside  is  positioned  outside  the  wind  tunnel.  The  wind  tunnel 
glass  reflects  the  light  emitted  by  the  pinholes  into  the  lens  of  the  camera. 
This  gives  the  appearance  of  two  dots  of  light  seemingly  suspended  in  the 
test  section  (Figure  2).  Fine  adjustment  of  the.  position  of  the  box  must  be 
made  so  that  both  dots  are  in  focus  and  do  not  interfere  with  the  dot  on  the 
nose  of  the  model. 

Measurements  taken  include  the  angle  of  the  camera  with  respect 
to  the  flow,  the  x,y,  and  z  coordinates  of  the  dots  with  respect  to  a  fixed 
point  (usually  the  position  of  the  nose  of  the  model  at  zero  angle  of  attack), 
and  the  distance  from  the  nose  of  the  model  to  the  jewel  bearing. 

With  the  camera  and  reference  dots  in  place,  the  model  is  given  a 
rolling  velocity  by  tangentially  directing  a  stream  of  compressed  air  on  its 
side.  The  model  is  held  in  place  by  means  of  a  rod  with  a  roller  bearing 
attached  to  the  end  which  fits  over  the  nose.  A  variable  frequency  strobe 
light  is  used  to  indicate  when  the  desired  initial  roll  rate  is  obtained.  The 
model  is  then  disturbed  and  allowed  to  pitch,  yaw,  and  roll  while  the  motion 
is  recorded  by  the  camera. 

Data  Reduction^ 

A  digital  computer  program  was  developed  to  convert  the  numerical 
information  provided  by  the  comparator  readings  of  the  test  films  into 

Euler  angles.  No  small  angle  assumptions  are  made;  therefore , the 
accuracy  of  the  test  data  is  increased. 

Considering  the  kinematic  relationships  of  the  test  set-up  with 
modified  Euler  angles  as  conventionally  defined  for  aeroballistic  axes,  the 
geometry  of  the  motion  of  the  model  is  shown  in  Figure  3. 

The  following  quantities  are  defined  as : 

R  =  distance  from  the  pivot  point  to  the  nose  of  the 
model 

yc,  zc  =  coordinates  of  nose  of  model  projected  on 
horizontal  plane 
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e  3.  Aeroballistic  Axis  System  and  Euler  Angles. 


,  z'Q  =  coordinates  of  nose  of  model  on  viewi 
of  camera 

y  =  inclination  of  line -of-sight  of  camera 
horizontal 
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the  air  flow  in  the  wind  tunnel,  and  the  Y,Z  plane  parallel  to  the  horizontal 
plane  of  the  wind  tunnel.  As  the  model  yaws  and  pitches,  the  aeroballistic 
axes,  x,y,z,  which  are  fixed  to  the  body  but  do  not  roll  with  it,  move 
through  the  yaw  angle ,  ip ,  and  the  pitch  angle ,  0  .  A  point  in  the  fixed  axis 
system  can  be  transformed  to  the  aeroballistic  axis  system  through  the 
matrix  equation 


(31) 


X 

'x' 

y' 

=  T 

Y 

z 

_Z 

where 


I 

!  COS0COS ip 

cos0sin0 

-sin0  \ 

T  = 

-sin  ip 

COS  ip 

0 

^  sin0cos^ 

sin0sin^ 

cos0  / 

(32) 


Therefore,  the  coordinates  of  the  nose  of  the  model,  which  are  (R,0,0)  in 
the  aeroballistic  axis  system,  are  expressed  by  Equations  (31)  and  (32)  as 


xc  =  RcosGcos^ 
yc  =  Rsin^cosQ 


z  =  -RsinQ 


(33) 


The  Euler  angles  can  then  be  found  from  equations  to  be: 


0  =  sin  ■*' 


ip  =  sin 


fr) 

\Rcos0/ 


(34) 


(35) 


Next,  the  relevant  camera  angle  relationships  are  derived.  As 
shown  in  Figure  4,  the  effect  of  the  camera  inclination  angle, y  ,  is  to 
rotate  the  Y,Z  plane  about  the  Y-axis  through  the  angle  y’  =  90-y  .  Since 
the  choice  of  where  the  viewing  plane  lies  is  arbitrary,  it  was  chosen  to 
coincide  with  the  plane  which  passes  through  the  origin  and  the  Y-axis  and 
is  inclined  at  the  angle  y'  to  the  Z-axis. 
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Y,Y' 

Figure  4.  Angle  of  Camera  Relative  to  Inertial  Axis  System 
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Thus,  a  point  in  the  fixed  axis  system  is  transformed  to  the  tipped 
axis  system  of  the  camera  through  the  matrix  equation 


x" 

X 

Y* 

ii 

Y 

Z\ 

_Z_ 

where  the  orthogonal  transformation 


!  siny 

0 

'■-cosy  A 

T  =. 

0 

1 

o  ■ 

Y  '  ! 

^cosy 

0 

sinyy 

(36) 


(37) 


Therefore,  the  coordinates  measured  from  the  film  frame,  y'c,  z'Q  can  be 
transformed  back  to  the  wind  tunnel  plane  through  the  inverse  of  Equation 


yc  =  yc 

zc  =  z’c/siny 


(38) 

(39) 


These  are  the  coordinates  necessary  to  determine' the  Euler  angles 
from  Equations  (34)  and  (35). 

It  is  necessary  to  derive  the  relationship  governing  the  data  on  the 
film  frame  and  unite  them  with  previous  relationships. 

Chances  are  that  the  plane  of  the  film  frame  in  thecamera  and  the 
viewing  plane  containing  the  reference  dots  do  not  coincide .  They  may  be 
assumed  parallel  due  to  the  largo  depth  of  focus  of  the  lens;  however,  a 
relative  translation  and/or  rotation  can  be  present  due  to  movement  of  the 
camera  during  a  test.  Consequently,  account  must  be  taken  to  eliminate 
this  possibility  of  error. 

If  the  axes  of  the  viewing  plane  are  defined  as  y' ,  z'.  and  the 


reference  points  are  R,  and  I^,  then 

Ri  *  (yj.  z'i> 

(40) 

r2  =  <y2-  z2> 

(41) 
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and  the  angle  of  rotation  is 


^  =  tan 


■1 


z2  -  Z>1 
y*o  -  y* 


(42) 


2  J  1 

If  the  axes  of  the  film  plane  are  defined  as  then 

R1  =  (yfl»zfl> 

R2  =  (yf2»zf2> 


(43) 

(44) 


Projection  of  the  film  plane  onto  the  comparator  yields  numerical  quantities 
for  Rp  R2,  and  the  nose  dot,  N,  in  terms  of  the  comparator  system 


R]_  =  (ni»n2) 


Ro  = 


(n3,n4) 


N  =  (n5>n6) 


(45) 

(46) 

(47) 


The  result  is  a  conversion  factor,  C,  from  comparator  units  of  measure¬ 
ment  to  the  actual  wind  tunnel  units: 


C  =  \ 


(nl  "  n3 )  +  (n2  "  H) 


D 


(48) 


where  D  is  the  measured  distance  between  the  reference  points. 

In  order  to  eliminate  comparator  bias  and  the  effect  of  a  translation 
by  the  camera  during  the  test,  the  comparator  origin  is  moved  to  the  left 
hand  (arbitrarily)  reference  point  by  letting 


R1 

N. 

O 

O 

II 

rH 

> 

•» 
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(49) 

R2 

*  <u2-v2>  •  <yfry(2-zf2'zfi) 

(50) 

N 

=  <u3'v3>  =  <yfl-y(3>  zf3  -  zf)> 

(51) 

The  angle  of  rotation  of  the  projection  of  the  film  with  respect  to  the 
comparator  can  now  be  computed  as 

€o  = 

A* 

&) 

(52) 
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Therefore,  the  net  rotation  to  be  considered  is 


(53) 

The  comparator  coordinates  can  now  be  transformed  to  the  wind 
tunnel  reference  system  by  a  rotation  through  the  angle  €  : 


/  cose 

sine\ 

u 

l-sine 

cose/ 

V 

where  u’  and  v*  are  the  desired  rotated  coordinates. 


(54) 


Finally,  these  coordinates  must  be  transformed  to  viewing  plane 
coordinates.  Since  the  position  of  the  origin  of  the  wind  tunnel  system  is 
the  position  of  the  model  nose  at  zero  angle  of  attack,  the  coordinates  of 
the  origin  with  respect  to  the  reference  points  can  be  computed  by  averag¬ 
ing  the  measured  relative  displacement  to  each  reference  point.  Letting 
Qy  and  Qz  be  these  coordinates,  then 

y'c  =  u3  ~  Qy 

z'c  =  »3  ■  Qz  <56> 

Substitution  of  these  quantities  into  Equations  (38)  and  (39)  and  the  results 
into  Equations  (34)  and  (35)  will  yield  the  desired  Euler  angles. 

Wind  Tunnel  Facility  and  Models' 

All  tests  were  conducted  in  contractor's  vertical  subsonic  wind 
tunnel  which  is  equipped  with  anti-turbulence  screens  and  a  16. 625- by 
16.625-inch  test  section.  The  maximum  wind  velocity  iS'  150  feet  per 
second;  however,  all  tests  were  run  at  approximately  45  feet  per  second. 
Employment  of  the  vertical  wind  tunnel  facilitated  the  use  of  a  low  friction, 
three -degrees -of- freedom,  jewel  bearing  model  mount.  The  model  is 
shown  mounted  in  the  wind  tunnel  in  Figure  5. 

Two  configurations  of  the  three  caliber  Army-Navy  Basic  Spinner 
Rocket  were  tested;  one  with  a  zero-degree  boattail  and  one  with  a  7.5- 
degree  boattail.  Both  models  were  2.5  inches  in  diameter,  constructed 
from  aluminum,  and  had  a  center  of  gravity  positioned  2.0  calibers  from 
the  nose.  The  7.5-degree  boattail  was  0.40  caliber  in  length.  Model 
exteriors  are  shown  in  Figures  6  and  7  and  other  physical  parameters  are 
listed  in  Table  I. 


Before  testing, the  models  were  balanced  in  still  air  to  eliminate  the 
effects  of  a  moment  due  to  their  own  weight.  An  adjustable  jewel  bearing 
mount  and  weight  ring  allowed  movement  of  both  the  center  of  gravity  and 
bearing  (Figure  8) . 

When  spinning  in  still  air  the  model  is  a  gyroscopic  pendulum.  A 
center  of  gravity  positioned  above  the  point  of  rotation  and  non-zero  angle 
01  attack  causes  a  precessional  motion  in  the  same  direction  as  the  rolling 
velocity.  Conversely,  a  center  of  gravity  positioned  below  the  point  of 
rotation  causes  precessional  motion  opposite  the  rolling  velocity.  If  the  two 
coincide,  the  model  is  neutrally  stable  and  no  precessional  motion  is 
caused.  Adjustment  of  the  weight  ring  and  bearing  position  was  made  until 


TABLE  I.  PHYSICAL  PARAMETERS  OF  TEST  MODELS 


the  last  condition  persisted.  The  jewel  bearing  was  positioned  exactly  the 
same  distance  from  the  nose  for  the  7.5-degree  boattail  model  as  it  was  for 
the  zero -degree  boattail  model.  Only  the  weight  ring  was  moved  in  balancing 
so  that  the  center  of  gravity  positions  of  the  two  models  were  identical. 


Friction  Analysis 

Frictional  effects  were  first  noticed  when  an  attempt  was  made  to 
coincide  the  center  of  gravity  of  the  model  with  the  position  of  the  jewel 
bearing.  Adjustments  had  been  made  so  that  no  angular  velocity  resulted 
over  a  time  span  of  15  seconds  (at  least  twice  the  duration  of  an  aerodynamic 
test);  however,  it  was  noticed  that  an  angular  displacement  decreased  in 
amplitude,  that  is,  when  displaced  the  model  tended  to  right  itself.  This 
way  theorized  to  be  a  frictional  effect  since  the  equations  of  motion  for  a 
gyroscopic  pendulum  predict  zero  damping.  (1)  Neglecting  this  phenomenon 
could  possibly  cause  serious  inaccuracies  in  the  dynamic  damping  factors. 
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Figure  8.  Model  Interior  and  Dimensions. 


Several  mathematical  models  were  conceived  to  handle  the 
problem;  however,  none  actually  gave  a  true  representation  of  what  was 
qualitatively  observed.  Consequently,  compensation  for  the  frictional 
effects  was  made  by  an  approximation.  Three -degrees-of- freedom  data 
was  taken  of  the  model  in  still  air.  The  data  was  then  reduced  to  obtain 
the  damping  factors  due  to  friction  (aerodynamic  forces  were  assumed 
negligible).  The  damping  factors  obtained  from  the  aerodynamic  tests 
were  then  corrected  through  the  following  equation: 


SECTION  V 


EXPERIMENTAL  RESULTS 


Data 

The  data  acquisition  technique  provided  two-axis  Euler  angles. 

Figure  9  shows  the  arms  of  tri-cyclic  motion  and  their  angular  rates.  No 
small  angle  assumptions  were  used  in  computing  the  9,  ft  angles  and  they, 
therefore,  contain  a  higher  degree  of  accuracy.  Samples  of  the  three - 
degrees -of-freedom  angular  data  are  shown  in  Figures  10  and  11  for  the 
zero-degree  and  7.5-degree  boattail  models,  respectively. 

! 

Inspection  of  Figure  10  shows  the  smaller,  faster  nutation  arm  rotat-  j 

ing  in  a  clockwise  direction  (the  same  direction  as  the  roll  rate).  A  large 
difference  can  be  seen  between  the  nutation  and  precession  frequencies. 

Comparison  of  Figures  10  and  11  shows  a  difference  in  the  relative 
sizes  of  the  nutation  and  precession  arms  of  the  two  models.  The  frequency 
ratio,  however,  is  approximately  the  same  for  both. 


Friction 

Reduction  of  the  data  films  taken  on  the  three -degrees -of-freedom 
motion  of  the  models  in  still  air  yielded  the  damping  characteristics  due  to 
friction  alone.  Results  indicate  that  the  precession  arm  was  strongly 
affected  while  negligible  effect  was  contributed  to  the  nutation  arm.  The 
still  air  damping  factors  were  small  but  of  the  same  order  of  magnitude  as 
the  dynamic  precession  damping  factors  for  both  models.  The  damping 
factors  for  the  nutation  arm  were  two  orders  of  magnitude  smaller  than  the 
dynamic  nutation  damping  factors;  hence,  they  were  neglected. 

The  value  of  Xpf  was  -0.023  for  the  zero-degree  boattail  and  -0.056 
for  the  7. 5 -degree  boattail.  Corrections  were  made  in  the  dynamic  pre¬ 
cession  damping  factors  using  these  quantities. 


Zero -Degree  Boattail  Analysis  I 

A  series  of  qualitative  wind  tunnel  tests  were  conducted  as  the  initial  { 

step  in  the  testing  procedure.  These  provided  an  estimate  of  the  aerodynamic  | 

characteristics  of  the  model  and  helped  establish  an  orderly  quantitative  ; 

testing  procedure.  \ 

! 

Initially  the  model  was  spun  up  to  approximately  900  rpm  (this  was  \ 
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close  to  the  maximum  roll  rate  attainable  before  the  model  would  separate 
itself  from  the  sting  due  to  vibrations).  The  model  was  then  disturbed  to 
excite  both  the  nutation  and  precession  modes.  The  resulting  motion  was 
comprised  of  a  low  frequency  precession  arm  having  an  amplitude  of  about 
nine  degrees  and  a  high  frequency  nutation  arm. 

The  precession  amplitude  increased  very  slightly  with  time,  indicat¬ 
ing  a  weak  instability.  On  the  other  hand,  the  nutation  amplitude  usually 
decreased  to  zero  before  the  completion  of  three  precession  cycles.  Also, 
at  900  rpm  it  was  virtually  impossible  to  disturb  the  model  in  such  a  way 
that  a  large  nutation  amplitude  could  be  attained.  Consequently,  it  was 
judged  that  the  nutation  arm  was  highly  stable  and  the  precession  arm  was 
slightly  unstable  at  this  roll  rate. 

Qualitative  observations  were  also  made  at  a  lower  roll  rate  to 
determine  possible  effects  of  roll  rate  on  the  dynamic  stability.  The  initial 
value  was  about  650  rpm.  No  lower  values  were  used;  this  insured  that  the 
model  possessed  adequate  gyroscopic  stability  throughout  the  test. 

At  650  rpm  an  initial  disturbance  resulted  in  a  larger  total  motion. 
The  initial  nutation  and  precession  amplitudes  were  greater  with  a  large 
change  in  the  size  of  the  nutation  arm.  The  nutation  amplitude  still  damped 
to  zero  but  much  more  slowly  than  at  the  higher  roll  rate. 

A  significant  change  was  noted  in  the  precession  mode  at  650  rpm. 
Insteady  of  the  slight  undamping  that  was  present  at  900  rpm,  the  amplitude 
slowly  decreased.  The  only  other  noticeable  difference  was  an  increased 
precession  frequency  and  a  decreased  nutation  frequency. 

The  performance  of  the  model  at  the  lower  roll  rate  and  higher  angle 
of  attack  indicated  the  presence  of  a  roll  rate  dependence  and  a  possible 
nonlinear  Magnus  moment. 

Since  the  purpose  of  the  test  program  was  to  determine  the  effects  of 
a  boattail  on  the  dynamic  stability  of  the  model,  an  attempt  to  determine 
roll  rate  effects  was  not  pursued.  However,  two  of  the  six  quantitative  tests 
of  the  zero-degree  boattail  were  conducted  at  a  lower  roll  rate  than  the  other 
four.  This  was  done  primarily  to  obtain  the  broadest  possible  range  of  com¬ 
plex  angles  of  attack  in  an  effort  to  sensibly  obtain  the  aerodynamic 
coefficients  as  nonlinear  functions  of  angle  of  attack. 

The  linear  quantitative  test  results  were  broken  up  into  two  sets  of 
three  each  -  the  thfee  lowest  and  the  three  highest  amplitude  funs.  All 
three  of  the  lower  amplitude  runs  [1(A),  1(B),  and  1(C)]  exhibited  a  stable 
nutation  mode  and  an  unstable  precession  mode.  The  higher  amplitude  runs 
[2(A),  2(B),  and  2(C)]  were  all  stable  in  both  modes.  The  stability  para¬ 
meters  are  shown  in  Figures  12  through  18  for  runs  1,  and  in  Figures  19 
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through  25  for  runs  2.  The  accuracy  of  the  stability  parameters  is  indi¬ 
cated  by  the  small  probable  errors -of- fit  plotted  versus  time  in  Figures  12 
and  19.  The  figures  also  include  the  nutation  and  precession  arms,  damp¬ 
ing  factors,  and  frequencies  plotted  versus  time. 
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Friction  caused  a  slight  roll  rate  decay  over  the  duration  of  a  test. 
The  effect  of  the  roll  rate  decay  was  time  varying  frequencies.  Because  of 
this  the  dynamic  damping  factors  were  corrected  on  the  basis  of  Equation 
(25). 

Inspection  of  the  numerical  values  for  the  dynamic  precession  damp¬ 
ing  factors  for  runs  2  showed  them  to  be  of  lesser  magnitude  than  the 
friction  damping  factors.  Consequently,  correcting  the  damping  factors  for 
friction  resulted  in  positive  values.  This  indicated  that  at  the  higher  angles 
of  attack  the  friction  was  overcoming  the  weak  precession  instability  .There¬ 
fore,  Figures  20  and  21,  which  show  decreasing  precession  arms  and 
positive  damping  factors  for  runs  2,  are  misleading.  The  precession  arms 
determined  by  the  fitting  routine  were  not  corrected  for  friction  before 
plotting,  whereas  the  damping  factors  were. 


Linear  aerodynamic  coefficients  were  calculated  from  Equations  (26),  1 

(27),  and  (28),  using  the  stability  parameters  obtained  from  the  "Wobble"  J 

program.  The  coefficients  are  plotted  versus  time  in  Figures  26  through  28  j 

for  runs  L  and  Figures  29  through  31  for  runs  2.  The  average  values  are  i 

listed  in  Table  II. 


TABLE  II.  LINEAR  COEFFICIENTS  (AVERAGE  VALUES) 


Run 

Cma  (1/Rad) 

Cmq  (1/Rad) 

cmp<*  O/Kad2) 

Ze  ro -degree 
Boattail 

1(A) 

1.0996 

-0.3472 

-0. 1795 

1(B) 

1. 1140 

-0.5913 

-0.1391 

1(C) 

1. 1051 

-0.8252 

-0. 1333 

2(A) 

1. 1299 

-1.1059 

-0.0957 

2(B) 

1. 1768 

-0.6422 

-0.0489 

2(C) 

1. 1262 

-0.2497 

-0.0440 

7. 5 -degree 
Boattail 

3(A) 

1. 1126 

1.2498 

0.2150 

3(B) 

1.0682 

1 . 3094 

0.2279 

3(C) 

1.0838 

1.4225 

0.0181 

3(D) 

1.  1194 

1.0500 

0.2492 

0.25 


28 


Figure  12.  Probable  Error-of-Fit  (Degrees)  versus  Mea 


Figure  13.  Precession  Arm  (Degrees)  versus  Mean  Time  (Seconds)  [  Runs  1(A),  1(B),  1(C)]  . 
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igure  14.  Dynamic  Precession  Damping  Factor  (One  Per  Second)  versus  Mean  Time  (Seconds) 
[Runs  1(A),  1(B),  1(C)]  . 


1(A) 


Mean  Time 

Precession  Frequency  (Radians/Second)  versus  Mean  Time  (Second)l  Runs  1(A),  1(B),  1(C)1  . 


-1(A) 

-1(B) 
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Figure  17.  Dynamic  Nutation  Damping  Factor  (One /Second)  versus  Mean  Time  (Seconds) 
t  Runs  1(A),  1(B),  1(C)]  . 


Figure  18.  Nutation  Frequency  (Radians/Second)  versus  Mean  Time  (SecondXRuns  1(A),  1(B),  1(C)]  • 


-  2(A) 

-  2(B) 

-  2(C) 


Figure  20.  Precession  Arm  (Degrees)  versus  Mean  Time  (Seconds) 
[Runs  2(A),  2(B),  2(C)]  . 
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Nutation  Arm  (Degree)  versus  Mean  Time  (Seconds)f  Runs  2(A),  2(B),  2(C)]  . 


e  25.  Nutation  Frequency  (Radians/Second)  versus  Mean  Time  (Seconds) 
[  Runs  2(A),  2(B),  2(C)  1  . 


Restoring  Moment  Coefficient  (One /Radian)  versus  Mean  Time  (Second) 
(Runs  L(A),  1(B),  1(C)]  . 


Figure  30.  Damping  Moment  Coefficient  (One /Radian)  versus  Mean  Time  (Second) 
[  Runs  2(A),  2(B),  2(C)]  . 


Of  primary  importance  is  the  near  constant  value  of  the  restoring 
moment  coefficient  for  any  particular  run.  One  of  two  possible  conclusions 
can  be  drawn:  either  the  angle  of  attack  range  is  small  so  that  the  nonlinear 
effect  is  not  seen,  or,  the  nonlinearity  itself  is  small  enough  to  be  neglected 
In  either  case,  the  constant  Cm  a  allows  the  roll  rate  to  be  calculated  as  a 
function  of  time  from  the  nutation  and  precession  frequencies  and  the  mo¬ 
ments  of  inertia  (Equation  5).  This  parameter  is  needed  for  the  nonlinear 
analysis. 

The  linear  Magnus  moment  coefficients,  Cm  ,  are  all  less  than 
zero.  From  the  Linear  Theory,  a  negative  CmPQ!  ™ul  have  a  stabilizing 

effect  on  the  nutation  arm  and  a  destabilizing  effect  on  the  precession  arm. 

The  damping  moment  coefficients,  Cmq  +0^^  »  are  all  also  less 
than  zero.  The  effects  of  a  negative  sign  are  to  staoilize  the  nutation  arm 
and  slightly  destabilize  the  precession  arm.  Combining  the  effects  of  the 
damping  moment  coefficients  and  Magnus  moment  coefficients  predict  a 
highly  stable  nutation  arm  and  an  unstable  precession  arm.  The  test  results 
support  this. 

Cmpa  and  Cmq  +  Cm<^  were  not  as  constant  over  a  given  run  as 

Cm  a  •  Indications  are  that  these  two  coefficients  are  much  more  nonlinear 
than  CmQf. 

Nonlinear  coefficients  were  determined,  using  the  procedure  outlined 
in  Section  III.  All  of  the  stability  parameters  from  the  six  test  runs  were 
evaluated  together  to  include  the  entire  angle  of  attack  range  covered  by  the 
tests,  thereby  giving  a  better  representation  of  the  nonlinear  characteristics 
of  the  coefficients.  The  nonlinear  coefficients  are  plotted  versus  angle  of 
attack  in  Figures  32  through  34  and  listed  in  Table  HI.  The  angle  of  attack 
range  and  corresponding  range  of  linear  coefficients  for  each  run  are  super¬ 
imposed  on  these  plots  to  show  how  well  the  linear  and  nonlinear  results 
compare. 

The  nonlinear  coefficient  results  are  in  good  agreement  with  the 
qualitative  observations.  Both  the  damping  and  restoring  moment  coefficients 
exhibited  hard  spring  characteristics  while  the  Magnus  moment  coefficient 
was  a  soft  spring  nonlinearity. 

As  expected,  the  nonlinear  term  in  Cma  was  small.  This  verifies 
the  prior  assumption  which  allowed  the  roll  rate  to  be  calculated.  Virtually 
no  effect  was  felt  by  the  nonlinearity  over  the  angle  of  attack  range  tested. 

^mq  +^m<i  was  3uite  nonlinear  and,  being  less  than  zero,  contri¬ 
buted  greatly  to  the  stability  of  the  nutation  arm  and  to  the  instability  of  the 
precession  arm.  These  effects  increase  with  angle  of  attack. 


48 


Figure  33.  Nonlinear  Damping  Moment  Coefficient  versus  Complex  Angle 
of  Attack  (Degrees)  (Runs  1  and  2). 


Nonlinear  Magnus  Moment  Coefficient  versus  Complex  Angle 
of  Attack  (Degree)  (Runs  1  and  2)  . 


k  r 

On  the  other  hand,  the  nonlinear  Magnus  moment  coefficient  has  an 
opposite  effect.  As  the  angle  of  attack  increases,  the  value  of  the  coefficient 
will  become  positive,  having  a  stabilizing  effect  on  the  precession  arm  and 
a  destabilizing  effect  on  the  nutation  arm.  These  effects  were  observed  in 
the  qualitative  tests. 

Due  to  the  nonlinear  effects  of  and  Cmq  +  ,  it  is  possible 

for  the  model  to  possess  dynamic  stability  at  higher  angles  of  attack.  How-  * 

ever,  a  critical  balance  must  be  attained  between  the  two  moments  along 
with  a  suitable  roll  rate.  The  roll  rate  would  have  to  be  high  to  reduce  the 
effects  of  Cmq  +  Cm(^  on  the  precession  arm. 

7. 5 -Degree  Boattail  Analysis 

Qualitative  tests  conducted  on  the  7. 5-degree  boattail  model  disclosed 
that  the  aerodynamic  characteristics  for  that  model  were  quite  different 
from  those  of  the  zero -degree  boattail  model. 

For  comparison  sake,  the  initial  roll  rate  was  about  900  rpm.  Low 
angle  excitation  of  the  model  resulted  in  a  motion  comprised  of  a  small  pre¬ 
cession  arm  of  decreasing  amplitude  and  a  small  nutation  arm  of  increasing 
amplitude. 

The  nutation  arm  undamped  rapidly  to  a  point  where  the  motion  was 
so  violent  that  the  model  would  separate  itself  from  the  sting.  In  some 
cases,  the  nutation  amplitude  would  build  up  from  near  zero  to  a  value 
appearing  to  be  twice  as  large  as  the  precession  amplitude  (approximately 
2.5  degrees).  In  some  cases  this  would  occur  before  the  completion  of  two 
precession  cycles. 

A  large  angle  disturbance  resulted  in  stability  in  both  the  nutation  and 
precession  modes.  A  precession  amplitude  of  approximately  eleven  degrees 
damped  at  about  the  same  rate  as  at  a  low  angle  disturbance.  The  nutation 
mode  was  extremely  stable.  It  was  impossible  to  create  a  nutation  disturb¬ 
ance  that  would  last  for  any  more  than  one -third  of  a  precession  cycle. 

If  the  model  was  disturbed  at  a  high  angle,  the  nutation  mode  would 
damp  quickly  along  with  a  moderately  damped  precession  mode.  The  total 
motion  would  then  finally  reach  the  small  angle  range  where  the  nutation 
amplitude  would  begin  to  increase  wildly.  It  was  then  concluded  from  the 
Linear  Theory  that  for  the  angle  of  attack  range  observed,  the  nutation  insta¬ 
bility  stemmed  from  an  improper  balance  between  a  positive  damping 
moment  and  a  positive  Magnus  moment. 

Quantitative  tests  dealt  only  with  the  small  angle  of  attack  range  only 
because  a  suitable  nutation  mode  could  not  be  attained  at  the  higher  angles. 
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A  larger  nutation  disturbance  may  have  been  possible  at  a  slower  roil  rate, 
but  the  need  for  like  conditions  to  compare  the  two  models  prevented  study 
in  this  area.  Four  test  runs  1  3(A),  3(B),  3(C)  and  3(D)  ]  were  conducted 
on  the  7.5-degree  boattail  model,  and  the  stability  parameters  were  ob¬ 
tained  by  analyzing  the  data  with  the  "Wobble"  program.  The  stability  para¬ 
meters  are  plotted  as  a  function  of  time  in  Figures  35  through  41.  The 
dyn  imic  damping  factors  were  corrected  for  varying  frequencies  due  to  the 
decaying  roll  rate  and  then  for  friction  as  with  the  zero-degree  boattail 
model  before  plotting. 

Run  3(C)  was  the  only  run  in  which  the  friction  damping  factor  was 
greater  in  magnitude  than  the  dynamic  precession  damping  factor.  This  con¬ 
dition  only  persisted  through  about  half  the  duration  of  the  test.  The  total 
motion  was  smaller  for  run  3(B)  than  for  the  other  runs  and  may  have  been 
a  contributing  factor. 

Linear  aerodynamic  coefficients  were  calculated  for  the  7.5-degree 
boattail  model  from  the  stability  parameters  obtained  from  the  "Wobble" 
program.  The  coefficients  are  plotted  in  Figures  42  through  44  and  the 
average  values  are  listed  in  Table  II. 

As  with  the  zero-degree  boattail  model,  CmQ,  remained  fairly 
constant  with  time  for  a  given  run.  Therefore,  Equation  (5)  could  again  be 
used  to  compute  the  roll  rate  for  the  nonlinear  analysis. 

The  average  Magnus  moment  coefficients  are  all  greater  than  zero. 
From  the  Linear  Theory  ,  a  positive  Magnus  moment  will  have  a  destabilizing 
effect  on  the  nutation  arm  and  a  stabilizing  effect  on  the  precession  arm. 


The  damping  moment  coefficients  are  also  greater  than  zero.  Positive 
values  will  destabilize  the  nutation  arm  and  stabilize  the  precession  arm . 


Added  together ,  the  effects  of  C 
stable  nutation  arm  and  a  stable 


mpa  anc!("mq+C 
precession  arm. 


mi 


predict  a  highly  un- 


From  the  large  angle  qualitative  tests,  the  possibility  of  highly  non 
linear  Magnus  and  damping  moment  coefficients  exists. 


Nonlinear  coefficients  were  obtained  by  evaluating  together  the 
stability  parameters  from  all  the  test  runs;  however,  the  angle-of-attack 
range  was  limited  so  a  possibility  of  error  exists.  The  nonlinear  coefficients 
are  plotted  versus  angle  of  attack  in  Figures  45  through  47  and  are  listed  in 
Table  III.  A  comparison  between  the  nonlinear  and  linear  results  is  also 
shown  on  the  figures. 


The  nonlinear  results  support  the  qualitative  observations;  the  restor¬ 
ing  and  Magnus  moment  possessed  hard  spring  characteristics  while  the 
damping  moment  was  of  the  soft  spring  type. 
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-  3(A) 

-  3(B) 


-  3(A) 

-  3(B) 


Dynamic  Precession  Damping  Factor  (One/Second)  versus  Mean  Time  (Second) 


-3(A) 


Mean  Time 

Nutation  Arm  (Degree)  versus  Mean  Time  (Second) [Runs  3(A),  3(B),  3(C),  3(D)] 
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-  3(C) 

-  3(D) 


uns  3(A),  3(B),  3(C), 


Nutation  Frequency  (Radian/Second)  versus  Mean  Time  (Second)  [Runs  3(A),  3(B) 
3(C),  3(D)]. 
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Damping  Moment  Coefficient  (One /Radian)  versus  Mean  Time  (Second) 
[Runs  3(A),  3(B),  3(C),  3(D)]  . 


Moment  Coefficient  (One /Radian2)  versus  Mean  Time  (Seconds) 


Nonlinear  Restoring  Moment  Coefficient  Cone /Radian  versus 
Complex  Angle  of  Attack  (Degree)  (Runs  3). 


Nonlinear  Damping  Moment  Coefficient  (One /Rad)  versus  Complex  Angle  of  Attack 
(Degrees)  (Run  3). 


Nonlinear  Magnus  Moment  Coefficient  (One/Radianz)  versus  Complex 
Angle  of  Attack  (Degree)  (Runs  3). 
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The  variation  in  C  with  angle  of  attack  was  again  found  to  be 
small.  The  smaller  angle -of- attack  range  reduced  the  nonlinear  effects 


even  more . 


^mq +^md 


was  highly  nonlinear,  changing  sign  from  positive  to 

negative  at  around  8.5  degrees.  It  was  this  that  caused  the  nutation  arm  to 
become  stable  at  the  higher  angles  of  attack. 


The  Magnus  moment  was  also  very  nonlinear  and  a  stabilizer  to  the 
precession  arm;  however,  its  destabilizing  effects  on  the  nutation  arm  are 
overshadowed  by  the  strong  damping  moment.  At  small  angles,  the  Magnus 
and  damping  moments  both  contribute  to  the  instability  of  the  nutation  arm 
and  the  stability  of  the  precession  mode. 
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SECTION  VI 


CONCLUSIONS 


The  three -degree s-of-freedom  dynamic  wind  tunnel  tests  of  the  three- 
caliber  Army-Navy  basic  spinner  projectile  model  demonstrated  the  feasibility 
and  reliability  of  obtaining  nonlinear  aerodynamic  stability  coefficients  from 
vertical  wind  tunnel  testing.  As  a  result  of  these  tests  on  the  zero -degree 
boattail  model,  it  is  concluded  that: 

•  The  model  exhibits  a  precession  instability  at  low  angles 
of  attack  due  primarily  to  a  negative  Magnus  moment. 

•  Dynamic  stability  is  possible  at  high  angles  of  attack. 

•  Due  to  the  nonlinear  characteristics  of  the  Magnus  and 
damping  moments  the  model  has  a  stable  nutation  mode  at 
all  angles  of  attack  (assuming  gyroscopic  stability). 

•  The  effect  of  the  nonlinearity  exhibited  by  Cm  „  was 
small. 

For  the  7.5  -  degree  boattail  model  the  tests  indicated  that: 

•  The  model  has  a  nutation  instability  at  low  angles  of 
attack  due  to  positive  Magnus  and  damping  moments. 

•  The  model  is  stable  in  both  nutation  and  precession 
modes  at  higher  angles  of  attack  due  to  the  nonlinear 
characteristics  of  the  Magnus  and  damping  moments. 

•  The  effect  of  the  nonlinearity  exhibited  by  CmQf  was 
small. 

It  is,  therefore,  concluded  that  addition  of  a  7.5-degree  boattail  to  the 
three-caliber  Army-Navy  basic  spinner  (c.g.  positioned  two  calibers  from 
nose)  destabilizes  an  otherwise  stable  nutation  arm  and  stabilizes  an  other¬ 
wise  unstable  precession  arm  at  low  angles  of  attack.  At  higher  angles  of 
attack  little  effect  is  noted  on  the  stable  nutation  arm  but  the  precession 
mode  becomes  stable. 
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section  vn 
RECOMMENDATIONS 


Because  of  the  cpmplex  relations  existing  between  the  damping  and 
Magnus  moments  of  the  projectile ,  additional  dynamic  typnel  tests  should 
be  conducted  to:  ”  ‘  * 

t  Determine  tp§  fluid  flow  properties  affecting  the  damping 
and  Magnus  moments, 

•  Establish  the  damping  apd  Magnus  characteristics  for  a 
wide  range  of  boattail  angles  and  sizes. 

•  Correlate  the  measured  damping  and  Magnus  coefficients 
with  the  flow  field  properties. 
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